Self-trapping transition in nonlinear cubic lattices 
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We explore the fundamental question about the critical nonlinearity value needed to dynamically 
localize energy in discrete nonlinear cubic (Kerr) lattices. We focus on the effective frequency and 
participation ratio of the profile to determine the transition into localization, performing several 
numerical simulations in one-, two-, and three-dimensional lattices. A simple criterium is developed 
- for the case of an initially localized excitation - defining the transition region in parameter space 
("dynamical tongue") from a delocalized to a localized profile. A general analytical estimate of the 
critical nonlinearity value for which this transition occurs is obtained. 



PACS numbers: 42.25.Dd, 42.65.Tg, 72.15.Rn 

Discrete nonlinear systems constitute a useful testbed 
to explore many interesting questions and properties of 
diverse branches of physics [IHl]. Different dimensions 
and topologies are possible, allowing the description of 
many different physical phenomena. Localized structures 
are natural solutions in this kind of systems, and they 
exist when a judicious balance between diffraction and 
self-focussing effects is reached. At low nonlinearity lev- 
els, initial excitations tend to diffract (expand) across the 
lattice due to the excitation of several extended states; 
at larger nonlinearity values, diffraction is inhibited and 
initial excitations remain localized. Thus, a very funda- 
mental question related to the necessary conditions to 
observe a dynamical transition from an extended pat- 
tern to a localized one appears. In other words, how 
large should the nonlinearity be in order to dynamically 
localize the excitation? The first attempts to solve this 
question where focused on small systems of just few sites. 
In Ref. [5 Kenkre and Campbell studied the self-trapping 
transition on a nonlinear cubic dimer, finding an exact 
value for the critical nonlinearity 7c = 4. Then, in Ref. [6] 
Molina and Tsironis explored systems ranging from 2 to 
100 sites. The general conclusion was that as the num- 
ber of sites increases, the transition approaches the dimer 
value being, therefore, not an increasing function of the 
lattice size. In Ref. ^ Dunlap et. al. studied this 
transition for a nonlinear cubic impurity embedded in 
a linear lattice and found that the critical nonlinearity 
was always lower than the one for a full nonlinear lat- 
tice. Bernstein et al. [8 compared the Hamiltonian of 
a single-site excitation with the one of a completely ho- 
mogeneous extended state to find that, for ID lattices, 
again, 7c = 4. In Ref. [9 Johansson et al. found numer- 
ically a lower critical value for ID lattices (7c ~ 3.45) 
computed in very large systems and for very long dy- 
namical time scales. This work also explored the transi- 
tion for ID binary lattices showing a smaller value com- 
pared to the monoatomic case. These authors mentioned 
that, by considering simple energy-balance conditions [8], 



a rough estimate can be obtained for simple lattices, but 
not for the binary case. A first attempt to find a more 
general criterium was done in Ref. [10] by comparing 7c 
with the minimum bound-state energy of a nonlinear im- 
purity problem. They found that, for nonlinear lattices 
of different dimensions and topologies, there is a kind 
of universal ratio of ~ 1.3 between these two energies. 
Very recently, Kevrekidis et al. [11] developed an analyt- 
ical criterium in terms of Hamiltonian comparisons. This 
analysis gave a sufficient, but not necessary, condition of 
7c = 4 (ID) and 7c = 7.3 (2D). 

In the present work we study the problem of local- 
ization in discrete nonlinear cubic (Kerr) lattices of dif- 
ferent dimensions and topologies. We develop a simple 
criterium to find the necessary nonlinearity to localize en- 
ergy in any lattice and verify this with direct numerical 
simulations. By studying the effective frequency and par- 
ticipation ratio evolution, we are able to identify clearly 
the localization-delocalization transition. Finally, we ex- 
tend our findings to binary (ionic or diatomic) lattices. 

The propagation of waves in nonlinear cubic lattices 
is well described by a discrete nonlinear Schrodinger 
(DNLS) equation pi44l: 

djUjrj^ _ _ A I I *P / \ 

- = ^nUn + VAfiUfi + -i\Ufi\ Ufi , (1) 

where Ufi corresponds to the wave amplitude at site 
n in a D-dimensional lattice of N sites. corre- 
sponds to the energy (propagation constant) at the n- 
th site, and V Afiu^i defines the linear interactions (cou- 
plings) between nearest-neighbor sites, z corresponds to 
the dynamical coordinate [1-4 . Without loss of gen- 
erality, we fix y = 1 and keep the nonlinear coeffi- 
cient "7" in order to study the necessary amount of 
nonlinearity to transit from a delocalized to a local- 
ized profile. Model ([T]) possesses two conserved quan- 
tities, the Norm Q = I'^nP, and the Hamiltonian 
H = Enbn\un\^^V{utAnun^c.c.)^{j/2)\un\^]. 
Extended linear (7 = 0) stationary solutions of model 
(IT| exist inside the linear spectrum {{Xbottom^ Kop} foi" 
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constant e^). Along this work, we will focus on the fol- 
lowing lattices: ID {Xtop = 2), 2D-Honeycomb (2D-h, 
Xtop — 3), 2D-square (2D-s, \to^ = 4), 2D-triangular 
(2D-t, \top = 6), and 3D-square (3D, Xtop = 6). We 
will study the self-trapping transition of delta-like input 
excitations: Ufi{z = 0) = (^n,no5 where no corresponds 
to a central position in the lattice (i.e., a bulk excita- 
tion). Considering that Q and H are dynamical con- 
stants, this input condition implies fixed values Qo = 1 
and Hq = Cfi^ +7/2. (A variation of 7 is equivalent to 
changing Q, the key parameter when thinking on exper- 
imental realizations p!2HT8] ). We use the participation 
ratio, defined as = Q^/ |^ as an indicator of the 

degree of localization of a wave-packet (e.g., = 1 for a 
single-site excitation, and R = N for an equally excited 
array) . 

To study the dynamical evolution and its relation 
with the stationary solutions of the system, we will 
consider an approximation implemented in Ref. |19j . 
that considers every instantaneous profile as an station- 
ary one. We assume stationary solutions of the form 
Un{z) = iifj exp (zAg^), being Ag the effective or average 
frequency of the profile. By replacing this in ([T]), mul- 
tiplying by and summing in all the sites along the 
lattice, we get an expression similar to the Hamiltonian. 
We rewrite this expression in terms of the fundamental 
quantities and parameters of model ([T]), obtaining that 
AgQ = + jQ^/{2R). However, Q = Qo and = Ho, 
therefore we obtain a closed and simple expression for 
the effective frequency: 

Ae = ..„ + i(l + |). (2) 

The dynamics is well described for this expres- 
sion: before the self-trapping transition, the wave- 
packet diffracts/delocalizes over the whole lattice and 
Ri ) ~ iV > 1 Ae( ) « eB„ + 7/2 {z 

max 

de- 
pends on the particular dimension, size, and topology of 
the lattice). On the other hand, after the self-trapping 
transition, localization dominates and R{zmax) ^ 1, im- 
plying that Xe{zmax) ^ ^fio + 7- Therefore, for nonlin- 
ear cubic lattices, we define a region in parameter- space 
-dynamical tongue- inside which all the dynamics is con- 
tained: 7/2 ^ Ae - Cfio ^ 7- 

By differentiating ([2| with respect to z we get 
dXe/dz = -{^/2R^)dR/dz. Therefore, below the tran- 
sition we must observe (on average) an evolution with 
an increasing R (diffusion) and a decreasing Ag, if 7 > 
0. Figure [l] shows an example for a 2D-square lattice 
where we observe that, before the self-trapping transi- 
tion (black, blue, and green curves), the effective fre- 
quency decreases from 7 to ^ 7/2 and that R increases 
from R = 1 to R ^ N ^ 1. For small values of 7 
(black and blue curves), the wave-packet approaches the 
lattice boundary at Zmax^ where is a maximum (see 
m Fig. [ijb) that for z > Zmax ^ 7, R decreases due to 
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Figure 1. (Color online) Example for a 2D-square lattice: (a) 
Ae and (b) R versus z for 7 = 0, 3, 7.17, and 7.18 represented 
by black, blue, green, and red lines, respectively. Dynamical 
evolution (\un,m(z)\'^) examples for (c) 7 = 3, (d) 7 = 7.17, 
and (e) 7 = 7.18. A/" = 31 x 31 sites and e^^ = 0. 



the reflections at the boundaries). Since Ag is not a con- 
served quantity of model ([T]), its evolution will present 
some oscillations around an average value [not observed 
in the scale of Fig. [ija)], implying that different, linear 
and nonlinear, frequencies are being excited. When the 
transition into localization takes place (red curves) the 
final average frequency becomes larger than 7/2, being 
outside of the linear band (Ag > Xtop)- The profile starts 
to excite more nonlinear frequencies reducing the excita- 
tion of linear modes and, as a consequence, it becomes 
localized (a similar process occurs in disordered lattices 
when studying nonlinear delocalization transitions [19^). 
Therefore, for our 2D-s example, self-trapping is observed 
for 7 ^ 7c = 7.18. The oscillation in Ag is accompanied 
by an oscillation (of opposite sign) of the i?- value [see red 
curves in Figs.[TJa) and (b)]. FigsjTJc)-(e) show some dy- 
namical propagation examples for different values of 7. 

We interpret Ag as an average frequency of the profile 
that includes all the modes excited during the propa- 
gation. To validate this interpretation, we compute a 
longitudinal discrete Fourier transform and obtain a nor- 
malized spectral density ^(A,7). Figure [2|a) shows this 
computation for a ID lattice. We see how a positive non- 
linearity modifies the original distribution of frequencies, 
exciting more modes closer to Xtop = 2. For 7 > 3.6 
[see the horizontal line in Fig. [2|a)], an emerging non- 
linear mode is strongly excited (brightest peak) and self- 
trapping starts to occur. In addition, we calculate a mean 
frequency, defined as (A) = A ^(A, 7) and plot, in Fig- 
ure [2|b), a comparison between this quantity and Ag. We 
see how both curves deviate from the 7/2-line in almost 




Figure 2. (Color online) Example for a ID lattice of = 401: 
(a) Spectral density versus 7 and A. (b) Xe{zmax) (black) and 
(A) (red) versus 7, inside the dynamical tongue (shaded area). 



the same region (7 > 3.6) [see Fig.[2|b)-inset]. In general, 
there is a very good agreement between both quantities. 
Only above the threshold, there is a larger disagreement 
due to the fluctuations in Ag. For higher dimensional lat- 
tices, we found a better agreement because fluctuations 
are reduced. These results validate our assumption of Ag 
as an effective and average frequency, becoming useful 
when studying different dynamical processes [19]. 

In Fig. [3l we collect our numerical results for differ- 
ent cubic lattices. Fig. [sj^a) shows that, for each lat- 
tice, all the dynamics is contained inside the "dynamical 
tongue" (shaded area), in agreement with our analytical 
prediction. When nonlinearity is small, the effective fre- 
quency is ~ 7/2 because the wave-packet have diffracted 
3) and (d)]. When the self-trapping occurs [see 
Ag increases in the direction of 7. For most of 
the studied lattices, we observe that the transition im- 
plies a rather abrupt jump in frequency {dXe/dj ^ 0). 
For one-dimensional systems, the transition is softer due 
to the absence of power thresholds to create a localized 
state 23 [21] . The excited nonlinear modes emerge from 
the top of the band [see Fig. [2|a)] and the interaction 
with linear modes is stronger (in fact, for 7 ~ 3.6, the av- 
erage frequency is smaller than Xtop implying that there 
is a mixture of linear and nonlinear excited modes). For 
dimension D > 2 there is a threshold to excite nonlinear 
modes, causing a fast jump in Ag above the band edge 
(the interaction of linear and nonlinear modes is weaker 
and transition becomes sharper). 

Fig. [sj^b) summarizes the results for the participation 
ratio normalized to Rq = R{j = O^Zmax)^ which cor- 
responds roughly to the size of the diffracted profile in 
the linear regime. We see that the self-trapping tran- 
sition implies a fast change in the participation ratio 
{dR/dj <C 0) from values of the order of the lattice 
size to values of the order of few lattice sites. This fig- 
ure shows that is a very useful parameter to trace the 
self-trapping transition, indicating an effective number of 
excited sites. This figure also shows that the transition 
is always abrupt, and sharper for dimensions larger than 
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Figure 3. (Color online) (a) Xe{zmax) and (b) R{zmax)/Ro 
versus 7 for ID (black), 2D-h (blue), 2D-s (green), 2D-t (red), 
and 3D (orange). Horizontal lines in (a) represent Xtop- Ver- 
tical dotted lines in (b) denote jc- For all lattices e^^ = 0. 



Table 1. Numerical results of 7c and Rc 



Lattice 


ID 


2D-h 


2D-S 


2D-t 


3D 




3.8 


5.46 


7.18 


10.8 


9.5 


Rc 


19 


10.1 


8.8 


9 


3.8 



one. By collecting the data from Fig. [3) we obtain the ap- 
proximate numerical 7c-values shown in Table [Tj We de- 
fine each one as the value of 7 for which Xtop = ^ei^max)- 
These values were obtained for fixed lattice sizes and 
fixed propagation distances (for different parameters, we 
obtain essentially the same 7c-values). Our aim is not to 
obtain an exact numerical value, because we believe that 
the relevant result is to obtain a numerical estimation 
for this transition (we should keep in mind that these 
transitions may occur and be observed in real physical 
systems, where an estimation range is more relevant [ij- 
4 ). In fact, if we use another quantity as an indicator 
(for instance, the space- averaged fraction of power [lOj), 
we observe the occurrence of the self-trapping transition 
in similar regions. 

Now, we tackle the complicated problem of getting 
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an analytical estimation of the critical nonlinearity (7c), 
for which the self-trapping transition occurs. In gen- 
eral, previous analysis observed that, at least, |7c| > 
\Kop — ^bottoml (which depends strongly on the particular 
dimension and topology of the lattice) . This criterium es- 
sentially says that when the nonlinear contribution (7Q) 
is larger than the linear one (approximately, the size of 
the band), the localization tendency will be more impor- 
tant than the diffr active one, and the wave-packet will 
tend to localize. This criterium has shown to be a very 
rough estimate of the actual value of 7c [9l [11] , which is 
indeed always an smaller value. Figures [l] and |3] show 
that the transition implies an abrupt change of the size 
of the wave-packet (proportional to R) and that the final 
effective frequency X{zmax) moves out of the band. Our 
hypothesis is quite simple: If the effective frequency is 
inside of the band, the profile is interacting with more 
linear modes (which are extended) and tend to delocal- 
ize. However, if the effective frequency is outside of the 
linear band, the wave-packet is able to excite a set of non- 
linear frequencies and only few linear modes, and tends 
to localize. In general, the transition will depend on the 
particular linear properties of the lattice but, also, on 
the sign of the nonlinearity. In this work, we concentrate 
only in the case 7 > 0, and on lattices possessing sym- 
metric linear bands (i.e., \Xtop\ = \XboUom\)- In this case, 
the 1 7c I -value is independent on the sign of the nonlin- 
earity and the data from Table [T] also applies for 7 < 0. 
For 7 > 0, we find that the critical nonlinearity is ob- 
tained when the effective frequency ^ coincides with the 
border of the linear band; i.e., when 

Ae = Xtop ^ 7c = 2 {Xtop - Cno) (^f^^^ ' 

In other words, we will be able to excite a nonlinear lo- 
calized state when the average frequency stays outside 
of the band and resonances with linear modes are re- 
duced or cancelled. We see that equation (|3| predicts 
a lower value of the critical nonlinearity than previous 
estimates. However, we get an extra parameter defined 
as Rc. We conjecture that this parameter corresponds to 
a critical size of the wave-packet for which self-trapping 
starts to occur, a kind of minimum volume for a pro- 
file to be considered as localized. In Table [l[ we present 
our computation of Rc for all the studied lattices using 
Eq. (|3| and the numerical 7c-values. Results show that 
this critical size decreases when the dimension increases, 
what is in agreement with the localization tendency of 
localized states in nonlinear cubic lattices. Close to the 
band, ID nonlinear stationary solutions are very broad, 
while for 2D lattices, and even more for 3D ones, non- 
linear stationary solutions are very localized above the 
norm threshold. For 2D lattices, we numerically found 
that this value is of the order of 10, leading to an estima- 
tion for 7c which is ~ 90% of the previous estimations. It 
is interesting to notice that Rc is nearly constant for two- 
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Figure 4. (Color online) Results for a binary 2D-square 
lattice, (al) [(bl)] shows R versus {7, Ae} for en^ — Ae 
[cno = 0]- (a2) [(b2)] shows \e{zmax) - versus {7, Ae} 
for efto Ae [enQ =0]. = 31 x 31 sites. 

dimensional discrete nonlinear cubic lattices, which are 
the most explored lattices nowadays in different contexts 
of experimental physics. 

Now we want to go further and prove the validity 
and robustness of our findings in more complex settings. 
For example, we study the critical values for binary lat- 
tices. We focus on a 2D-square array with site-energies 
C2n+i = and 62^^ = Ae, being Ae the site-energy con- 
trast. In this case, the position of the input excitation 
plays a fundamental role in the determination of the tran- 
sition. The factor "Atop — e^^" in expression (|3| implies 
a lower critical value for a larger e^^ . After a straightfor- 
ward computation, we obtain the band edge frequency 
Xtop = 0.5(Ae + V64 + Ae2). Taking R^ = 10, we get the 
following estimate for the critical nonlinearity 

7c ^ 0.9 (V64 + Ae2 + Ae - 2e^,) . (4) 

Therefore, for Cfi^ = Ae we expect a reduction of 7c if 
the contrast increases (Ae 1 ^ 7c ^ 0). On the other 
hand, for e^^ = we expect an increment of the critical 
nonlinearity if the contrast increases. Figure [4] collects 
our findings for this case. The density plots (al) and 
(bl) show the participation ratio versus 7 and Ae; dark 
(bright) color means an smaller (larger) i?-value, and the 
estimate Q is plotted by orange straight lines. We see an 
almost perfect agreement between the numerics and the 
analytical estimation, what validates the used Rc-yslue. 
Figures (a2) and (b2) show the transition in terms of the 
effective frequency. We clearly see the opposite tenden- 
cies depending on the input site-energy. Xe{zmax) — ^Hq is 
effectively 7/2 below the transition (bright fiat surface), 
and it increases in the direction of 7 above 7c. The tran- 
sition in participation ratio and frequency is very abrupt 
and perfectly described by Q. 

In conclusion, we have studied the problem of the self- 
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trapping transition for different nonlinear cubic lattices. 
We showed, theoretically and numerically that, by con- 
sidering a single-site input excitation, all the dynamics is 
contained in a very precise parameter-region in the shape 
of a "dynamical tongue". We have found approximate 
numerical values for the critical nonlinearity where this 
transition occurs. We showed that the effective frequency 
can be understood as an average quantity, giving a good 
insight of the frequencies participating in the dynamics. 
We analytically developed a new criterium to predict self- 
trapping, that makes use of an extra parameter (effective 
volume) that depends strongly on the dimension of the 
system. Binary lattices are well described by our analyt- 
ical prediction. 
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